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Planetary ephemerides are a very powerful tool to constrain deviations from the theory of General 
Relativity using orbital dynamics. The effective field theory framework called the Standard-Model 
Extension (SME) has been developed in order to systematically parametrize hypothetical violations 
of Lorentz symmetry (in the Standard Model and in the gravitational sector). In this communication, 
we use the latest determinations of the supplementary advances of the perihelia and of the nodes 
obtained by planetary ephemerides analysis to constrain SME coefficients from the pure gravity 
sector and also from gravity-matter couplings. Our results do not show any deviation from GR 
and they improve current constraints. Moreover, combinations with existing constraints from Lunar 
Laser Ranging and from atom interferometry gravimetry allow us to disentangle contributions from 
the pure gravity sector from the gravity-matter couplings. 

PACS numbers: 04.50.Kd,04.80.Cc,11.30.Cp 


I. INTRODUCTION 

The Solar System has proven to be an efficient labora¬ 
tory to discover new phenomena from gravitational ob¬ 
servations. Historically, one can mention the discovery 
of “dark” components (such as the planet Neptune pre¬ 
dicted by Le Verrier) or evidence towards non-Newtonian 
gravity theories (for example the perihelion advance of 
Mercury which pointed towards General Relativity - 
GR). The Solar System remains the most precise lab¬ 
oratory to test the theory of gravity, that is to say GR. 

Gonstraints on deviations from GR can only be 
obtained in an extended theoretical framework that 
parametrizes such deviations. The constraints that are 
obtained from observations are framework-dependent. In 
the past decades, two frameworks were widely used in the 
literature at the scale of the Solar System, namely the 
Parametrized Post-Newtonian (PPN) formalism [1] and 
the fifth force framework [2]. Stringent constraints have 
been obtained for these formalisms [1, 3-8]. More re¬ 
cently, other phenomenological frameworks have been de¬ 
veloped like the Standard-Model Extension (SME). The 
SME is an extensive formalism that allows a systematic 
description of Lorentz symmetry violations in all sec¬ 
tors of physics, including gravity [9-11]. Violations of 
Lorentz symmetry are possible in a number of scenarios 
described in the literature. While some early motivation 
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came from string theory [12], Lorentz violations can also 
appear in loop quantum gravity, noncommutative field 
theory and others [13, 14]. The SME is an effective field 
theory aiming at making phenomenological connections 
between fundamental theories and experiments. 

In particular, a hypothetical Lorentz violation in the 
gravitational sector naturally leads to an expansion at 
the level of the action [11, 15] which in the minimal SME 
writes 

5grav = J {R-UR + s^^'Rl, + 

( 1 ) 

with G the gravitational constant, g the determinant of 
the metric, R the Ricci scalar, the trace-free Ricci 
tensor, Cap^iv the Weyl tensor and u, and 
the Lorentz violating helds. To avoid conflicts with the 
underlying Riemann geometry, we assume spontaneous 
symmetry breaking so that the Lorentz violating coef¬ 
ficients need to be considered as dynamical fields. The 
last part of the action S' contains the dynamical terms 
governing the evolution of the SME coefficients. In the 
linearized gravity limit, the metric depends only on u and 
which are the vacuum expectation value of u and 
s'"' [15]. The coefficient u is unobservable since it can 
be absorbed in a rescaling of the gravitational constant. 
The so obtained post-Newtonian metric differs from the 
one introduced in the PPN formalism [15]. In addition 
to the minimal SME action given by Eq. (1), there exist 
some higher order Lorentz-violating curvature couplings 
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in the gravity sector (non-minimal SME) [16] that have 
been constrained by short range experiments [17]. These 
terms are not considered in this communication. 

In addition to Lorentz symmetry violations in the pure- 
gravity sector, violations of Lorentz symmetry can also 
arise from gravity-matter couplings. In [18], it has been 
shown that gravity-matter couplings violation of Lorentz 
symmetry can be parametrized by the following classical 
point mass action 

‘^mat — dA ^ d” 'll (Ueff) 5 ( 2 ) 

where is the four-velocity of the particle, m is its mass 
and and (oeff)^ are Lorentz violating fields. In this 
action, spin-coupled Lorentz violation is effectively set 
to zero. The new fields c^'' and (oeff)^ depend on the 
composition of the point particle [18]. This modification 
of the action produces two different types of effects: (i) 
a modification of the way gravity is sourced and (ii) a 
violation of the three facets of the Einstein Equivalence 
Principle. The first effect will result in a modification 
of the space-time metric solution of the field equations. 
Modifications of the metric in the linearized approxima¬ 
tion depend on (afg) coefficients, the background values 
of the coefficients (oeff)^ from the source body [18]. On 
the other hand, the violation of the equivalence princi¬ 
ple generated by the action (2) leads to a deviation from 
the geodesic motion depending at first order on the co¬ 
efficients and (a^)^, the background values of the 
Lorentz violating fields of the test mass. 

Up to now, several studies have constrained the pure- 
gravity SME coefficients like for example Lunar Laser 
Ranging [19], atom interferometry gravimetry [20, 21], 
short range experiment [22], planetary orbital dynam¬ 
ics [23], Gravity Probe B [24] and recently binary pul¬ 
sars [25]. The coefficients are currently poorly 

constrained by [26-29]. On the opposite, some of the 
coefficients are severely constrained (see for exam¬ 
ple [26, 27, 30, 31]). A list of current constraints on all 
SME coefficients can be found in [32]. In this study, we 
will concentrate on the impact of and coeffi¬ 

cients on planetary orbital dynamics and neglect the 
coefficients and leave them for future work. 

In this communication, we show that planetary or¬ 
bital dynamics can be used to derive stringent constraints 
on the SME coefficients. Indeed, SME modifications of 
gravity induce a secular variation of some orbital el¬ 
ements [15, 18] such as the longitude of the ascend¬ 
ing node and the argument of perihelia. These varia¬ 
tions are introduced in Sec. II. In Sec. Ill, we compare 
these variations with the present level of residuals com¬ 
ing from INPOPlOa (Integrateur Numerique Planetaire 
de rObservatoire de Paris) ephemerides [33]. We use a 
Bayesian inversion to infer the posterior probability den¬ 
sity function (pdf) on the SME coefficients. From the 
pdf, we estimate correlations between the coefficients. 
We estimate realistic confidence intervals and also de¬ 
termine linear combinations of the SME coefficients that 


can be determined independently from planetary orbital 
dynamics. In Sec. IV, we combine our results with pre¬ 
vious results obtained by Lunar Laser Ranging analysis 
and atom interferometry gravimetry. Finally, in Sec. V, 
we discuss our obtained results and present several ideas 
that may improve the current analysis. 


II. EFFECTS OF SME ON ORBITAL 
DYNAMICS 


In the linearized gravity limit, the gravity sector of 
the minimal SME is parametrized by a symmetric trace 
free tensor s^’' and by a scalar u that is unobservable 
since it corresponds to a rescaling of the gravitational 
constant [15]. Furthermore, the matter-gravity coupling 
is parametrized amongst others by the (oeg)^ coefficients 
which depend on the composition of the different bod¬ 
ies. The components of these coefficients depend on the 
observer coordinate system. The standard frame used in 
the SME formalism labeled by (T, X, Y, Z) is comoving 
with the Solar System, the spatial axes are defined by 
equatorial coordinates (see Fig. I of [15]) and the origin 
of time is given by the time when the Earth crosses the 
Sun-centered X-axis at the vernal equinox. The plan¬ 
etary orbital elements are defined with respect to the 
ecliptic coordinate system. The two coordinate systems 
differ by a rotation TZ of angle e = 23.44 ° (the Earth 
obliquity) around the X axis. Therefore, the transfor¬ 
mation of the tensor is given by sb = TViTVjS^'’ and 
gOi _ where capital letters refer to the equato¬ 

rial reference system and lower case letters refer to the 
ecliptic one. Similarly, the transformation of the (heff)^ 
vector is given by (a^g)* = Tdj {a^eY■ 

SME modifications of gravity induce different types 
of effects (for an extensive review, see [15, 18]). Two 
important effects can have implications on planetary 
ephemerides analysis: effects on the orbital dynamics and 
effects on the light propagation. Simulations using the 
Time Transfer Formalism [34] based on the software pre¬ 
sented in [35] have shown that only the and (hgg) 
coefficients produce a non-negligible effect on the light 
propagation (while it has impact only at the next post- 
Newtonian level on the orbital dynamics [15, 18]). Since 
in this analysis we concentrate on orbital dynamics, these 
coefficients are not considered and will be neglected. This 
can safely be done since the signatures from the s’^'^ and 
(hgg)^ coefficients on the light propagation are similar 
to the logarithmic standard Shapiro delay, which is not 
correlated to orbital dynamics effects. 

The equations of motion in the SME formalism are 
given in [15, 18]. Neglecting the contributions, the 
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two-body equation of motion reads 
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where Gjv is the observed Newton constant, M = mi + 
m 2 is the total mass of the two bodies, 5m = m 2 — mi 
is the difference of the two masses, P = r{ — r 2 is the 
relative position of the two masses and 

n- = , (4) 


with Ni 2 the number of particles of species w in the body 
1,2. The coefficient is completely unobservable in 
this context since absorbed in a rescaling of the gravita¬ 
tional constant (see the discussion in [15, 24]). The coef¬ 
ficient can also be absorbed in a rescaling of the 

gravitational constant that depends on the composition 
of each planet [18]. In this context, one would observe 
a different Gn with the different planets. Nevertheless, 
this effect is expected to be very small [18] and would not 
produce any supplementary advances of the perihelia and 
of the nodes and therefore is neglected in this analysis. 

In Eq. (3), the sums on w need to be done on the 
electrons, protons and neutrons. In the case of a Sun- 
planet system, we have M = m^ + mQ « toq, 5m = 
mQ — mp « mQ and 713 = — Nq « —Nq. The fact 

that we are neglecting means that we are neglect¬ 
ing effects produced by the violation of the universality 
of free fall. Under these assumptions, the equations of 
motion depend on 


N, 


w 

« - 0.9a {dlg^y - 0.1a , 


(5a) 

(5b) 


where we used a simple model for the composition 
of the Sun characterized by iVQ/mQ = N^/mQ « 

0.9(GeV/c^) ^ and N^/mQ « 0.1(GeV/c^) ^ as de¬ 
scribed in [18] (with c the speed of light in vacuum). In 
this paper, a (a“g)'^ is always expressed in GeV/c^ and 


(aS-)' = &)’ + (5?,)' . 


( 6 ) 


Using the Gauss equations, secular perturbations in¬ 
duced by SME on the orbital elements can be computed 
similarly to what is done in [15, 23]. The two orbital el¬ 
ements needed for our analysis are the longitude of the 
ascending node U and the argument of the perihelion w. 


The secular change in these two elements is given by 
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(7b) 


where a is the semimajor axis, e the eccentricity, i the 
orbit inclination (with respect to the ecliptic), n = 
{GN'm.Q/P)^^^ is the mean motion and e = 1 —( 1 — 

In all these expressions, the coefficients for Lorentz vio¬ 
lation with subscripts P, Q, and k are understood to 
be appropriate projections of along the unit vectors 
P, Q, and k, respectively. For example, Sq = PSq, 
Spp = The unit vectors P, Q and k define the 

orbital plane 


P = (cos U cos UJ — cos i sin U sin uj) (8a) 

-I- (sin U cos UJ + cos i cos U sin uj) Cy + sin i sin we),, 
Q = — (cos U sin w -I- cos i sin U cos uj) Sx (8b) 

-I- (cos i cos fl cos UJ — sin U sin uj) Cy + sin i cos ujCz , 
k = sin i sin Ueq, — sin i cos Ocy -I- cos ie^ , (8c) 

where Cx,y,z define the basis of the ecliptic reference sys¬ 
tem. The relations (7) are generalizations of Eqs. (168- 
171) from [15] that do not include the (d'^g)^ terms. 


III. ANALYSIS AND RESULTS 

Planetary ephemerides analysis uses an impressive 
number of different observations to produce high accu¬ 
rate planetary and asteroid trajectories. The observa¬ 
tions used to produce ephemerides comprise radioscience 
observations of spacecraft that orbited around Mercury, 
Venus, Mars and Saturn, flyby tracking of spacecraft 
close to Mercury, Jupiter, Uranus and Neptune and op¬ 
tical observations of all planets [4, 6 , 7, 33, 36-43]. Esti¬ 
mations of supplementary advances of perihelia with the 
Russian Ephemerides of Planets and the Moon (EPM) 
are presented in [ 6 , 42]. The INPOP ephemerides have 
produced estimations of supplementary advances of per¬ 
ihelia and nodes. Tab. I gives estimations obtained by 
INPOPlOa [33] on supplementary longitude of nodes ft 
and on supplementary argument of perihelia^ uj. 


^ In [.jS], LJ is noted uj which is commonly used for the longitude 
of the perihelion but the estimated values correspond to supple¬ 
mentary argument of perihelia and not to longitude of perihelia 
(usually noted by zu) [44]. 
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TABLE 1. Values of supplementary longitude of nodes and 
argument of perihelia estimated by INPOPlOa (see Tab. 5 
from [33]). These values are estimated in [33] as the interval 
in which the differences of postfit residuals are below 5 %. 

Planet 

Q (mas X cy ^) 

oj (mas X cy ^) 

Mercury 

1.4 ± 1.8 

0.4 ± 0.6 

Venus 

0.2 ± 1.5 

0.2 ± 1.5 

EMB 

0.0 ±0.9 

-0.2 ±0.9 

Mars 

-0.05 ±0.13 

-0.04 ±0.15 

Jupiter 

-40 ± 42 

-41 ± 42 

Saturn 

-0.1 ±0.4 

0.15 ±0.65 


Since and (a^)^ do not play any role in the orbital 
dynamics and s^’' is trace free, the observations depend 
on 8 independent fundamental coefficients: 

sQ = sXX^^YY_2gZZ^ ^XY^ ^YZ gTJ (^^ese 

coefficients will be denoted as pi in the following). In 
this communication, we perform a Bayesian inversion to 
infer knowledge on these 8 independent coefficients us¬ 
ing a Monte Carlo Markov Chain (MCMC) algorithm. 
The approach is very similar to the one used for binary 
pulsar data [25]. The observations are assumed to be 
independent and the errors to be normally distributed. 
The pdf describing the likelihood (i.e. the probability to 
obtain observations Oi given certain values of the SME 
coefficients pk) is given by 


L(Oi|pi,P 2 , • ■-Pn) = cste (9) 


where the is computed by 


X^ = 


E 

pi 


(Wpi.SME(p/c) — d>p/jNPOp)^ 
(ppl,SMEiPk) — INPOP j 




( 10 ) 


where the index pi of the sum is running over the six dif¬ 
ferent planets from Tab. I, fIpjjNPOP) WpzyNPOP and the 
corresponding a are from Tab. I and where d;p;_sME(Pfc) 
and fIp; 3 ME(Pfc) are simulated values depending on the 
SME coefficients by (7). The posterior pdf of the SME 
coefficients is given by 


P{pi,P2, . . .Pn|0*) = C L{Oi\pi, . . .p„)7r(pi, ...Pn), 

( 11 ) 

where 7r(pi,.. .p„) = 7r(pi)... 7r(p„) is the prior pdf on 
the SME coefhcients pk and C a constant. We use a uni¬ 
form prior pdf on the SME coefficients and the MCMC 
algorithm used is a standard Metropolis-Hasting algo¬ 
rithm [45]. We run the Metropolis-Hastings sampler un¬ 
til 10® samples have been generated. The convergence 
of the MC is ascertained by monitoring the estimated 
Bayesian confidence intervals of the parameters. Finally, 
to diminish the effect of the starting configuration, we 
discard the first 1000 samples. 


The marginal pdf of a single SME coefficient pj is given 

by 

PiPj\Ot) = J dpi J dp2 ■ • ■H(pi,... ,p„10i), (12) 

where the integrals are performed over all the SME coef¬ 
ficients Pk except Pj. 

A first run shows that the coefficients of our model 
are highly correlated, see Fig. 1. We have used the cor¬ 
relation matrix estimator to assess the strength of the 
parameters correlations, see Tab. III. These correlations 
are mainly due to the fact that all planets have very sim¬ 
ilar, low inclination, orbital planes. Nevertheless, we can 
produce marginal ID posterior distribution for each of 
the 8 SME coefficients. The histograms corresponding to 
these distributions are presented in Fig. 1. The corre¬ 
sponding Bayesian confidence intervals are presented in 
Tab. II. 


TABLE II. Estimations of the SME coefficients. These esti¬ 
mations are still correlated and the correlation matrix is given 
in Tab. III. The uncertainties correspond to the 68% Bayesian 
confidence levels of the marginal pdf. 


SME coefficients Estimation 


-XX -YY 

s — s 

(-0.8 ±2.0) X 10"^° 

-XX ^ -YY _ 

2 (-0.8 ±2.7) X 10"^° 


(-0.3 ± 1.1) X 10"^° 


(-1.0 ±3.5) X 10"“ 


(5.5 ±5.2) X 10"^^ 

qTX 

(-2.9 ±8.3) X 10"® 

CTY 

^0 

(0.3 ± 1.4) X 10"® 

qTZ 

_ 

(-0.2 ±5.0) X 10"® 


TABLE III. Estimations of the correlations coefficients be¬ 
tween the different SME coefficients: 


-rz qTX 
’ ’ ‘^0 

qTY 
, *^0 

and S'g^. 


1 




0.99 

1 



0.99 

0.99 

1 


0.98 

0.98 

0.99 1 


-0.32 

-0.24 ■ 

-0.26 -0.26 1 


0.99 

0.98 

0.98 0.98 -0.32 

1 

0.62 

0.67 

0.62 0.59 0.36 

0.60 1 

-0.83 

-0.86 - 

-0.83 -0.81 -0.14 

-0.82 -0.95 1 


Another approach (based on the first run) to avoid 
highly correlated coefficients is to find the independent 
linear combinations of the SME coefficients that can 
be determined by planetary ephemerides analysis. This 
can be done numerically by performing a normalized 
Cholesky decomposition of the covariance matrix 

C = K'^D^K, 


(13) 
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FIG. 1. 2D marginal posterior pdf (useful to assess the corre¬ 
lations). On the 2D plots, the blue dotted contours represent 
the 67 % Bayesian confidence area, the red continuous con¬ 
tour represent the 95 % Bayesian conhdence area and the 
dashed green contours represent the 99.7 % Bayesian confi¬ 
dence area. The histograms represent the marginal pdf of the 
SME coefficients. 


where C is the covariance matrix of the SME coefficients 
estimated from our first run, K is an upper triangular 
matrix whose diagonal elements are unity and is a 
diagonal matrix. Then the linear combinations b of the 
fundamental SME coefficients (noted p) given by 

b = K~^p , (14) 

with K~^ the inverse of the transpose of K, can be deter¬ 
mined completely independently by the analysis of plan¬ 
etary orbital dynamics. In our case, this Cholesky de¬ 
composition {K~'^) is given by 

b, = - s^y) , (15a) 

62 =-1.376i(15b) 

63 = -0.156i-0.31s'5-hs^^, (15c) 

64 = 0.0135i -h 0.064 s' 3 - 0.48s^^ -h , (15d) 

65 = 0.2661 - 0.31s' 3 O.Sls'^’^ - 1.67s^^ 

(15e) 

be = -35.561 -h 9.355*3 _ 22.675^"^ - 33.955^^^ 

-b7.83s^^ -f Sq , (15f) 

67 = 1641.46i - 2101.15*3 -h 4939.95^^ - 8846.85^^ 
-b4810.65^^ - 0.895^ -h , (15g) 

bs = 44.561 -h 47.15*3 - 580.15^'^ -f 1041.35^^ 

-b231.55^^ -h 3.435^ -h 2.56^^ -h , (15h) 


with the expression of Sq given by Eq. (5) . We can now 
use the linear combinations 6 ^ as fundamental parameters 
for our analysis. Performing a new MC run (using the 
same prior and likelihood as previously), we show that 
they can be estimated without any correlation. This can 
be seen in Fig. 2 where the 2D marginal posterior pdf on 
the bi combinations are presented. More quantitatively, 
the computation of the correlation matrix shows that the 
bi combinations are completely decorrelated by planetary 
ephemerides analysis since the absolute values of the cor¬ 
relation parameters never exceed 0.03. The ID posterior 
pdf of the bi combinations are also represented in Fig. 2. 
The estimated mean and standard deviation are given 
in Tab. IV. The obtained uncertainties are much smaller 
than those given in Tab. II. 


TABLE IV. Estimations of the independent linear combina¬ 
tions bi of the SME coefficients. The expressions of the com¬ 
binations bi are given by Eqs. (15). The uncertainties corre¬ 
spond to the 68% Bayesian conhdence levels of the marginal 
pdf. 


SME linear combinations 

Estimation 

61 

(-0.8 ±2.0) X 10"^*’ 

^2 

(2.3 ±2.3) X 10"“ 

63 

(3.0 ±9.7) X 10"^2 

64 

(0.2 ± 1.1) X 10"^^ 

65 

(-0.3 ±2.4) X 10"^® 

be 

(0.2 ± 1.1) X 10"® 

br 

(-0.6 ±2.3) X 10"® 

bs 

(0.3 ± 1.7) X 10"® 


We want to emphasize the fact that the results from 
both approaches presented above are completely equiv¬ 
alent. They are two ways to represent the same results. 
One is free to choose which approach is more appropri¬ 
ate: to work with the fundamental SME coefficients de¬ 
termined by Tab. II at the price of including the covari¬ 
ance matrix (or equivalently the correlation matrix from 
Tab. Ill) in the analysis or to work with uncorrelated 
linear combinations of the SME coefficients that are de¬ 
termined by Tab. IV. The results provided by both ap¬ 
proaches describe the same physical information. There¬ 
fore, they are completely equivalent. 


IV. COMBINATION WITH LUNAR LASER 
RANGING AND ATOM INTERFEROMETRY 
GRAVIMETRY 


It is interesting to combine the results obtained in the 
last section with constraints available in the literature. In 
particular. Lunar Laser Ranging (LLR) data have been 
used to constrain the pure gravity sector of SME [19]. 
Similarly, atomic gravimetry data have also been used to 
constrain the s^‘' coefficients [20, 21]. We will first com¬ 
bine our results from Sec. Ill with LLR results to produce 
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ibi ^2 ^3 ^4 ^5 ^6 ^7 ^6 


FIG. 2. 2D marginal posterior pdf (useful to assess the corre¬ 
lations) of the linear combinations bi of the SME coefficients 
given by Eqs. (15). On the 2D plots, the blue dotted contours 
represent the 67 % Bayesian confidence area, the red continu¬ 
ous contour represent the 95 % Bayesian confidence area and 
the dashed green contours represent the 99.7 % Bayesian con¬ 
fidence area. The ID histograms represent the marginal pdf 
of the SME linear combinations bi. 


constraints on the SME pure gravity sector alone. This 
will highlight the improvement brought by the planetary 
ephemerides data. In a second step, we will consider both 
the pure gravity sector and the gravity-matter couplings 
coefficients. We will demonstrate that the combination 
of planetary ephemerides data, LLR data and atom in¬ 
terferometry gravimetry data allows one to completely 
disentangle all the SME coefficients and (a^g)'^. 

The procedure to combine different types of analysis is 
standard and consists of performing a global least squares 
fit of all the estimations available. Obviously, the plan¬ 
etary estimations given by Tab. II are not independent. 
To take into account the correlation between the coef¬ 
ficients estimated in Sec. Ill, we use the parameter co- 
variance matrix from Tab. Ill as a weight in the least 
squares fit. Similarly, the coefficients estimated in the 
LLR analysis are weighted by their standard deviation 
in the least squares fit. Since no covariance matrix can 
be found in the literature, we assume these estimations 
to be independent (this corresponds to a worst case sce¬ 
nario). Instead of working with results given in Tab. II 
that are correlated, we can equivalently use the linear 
combinations given by Eqs. (15) and we then use the es¬ 
timated standard deviations from Tab. IV to weight the 
least squares fit. In that approach, the weight matrix 
in the fit is diagonal. We insist on the fact that both 


approaches lead to the same results. In the following we 
provide the mean and the standard deviation of the SME 
coefficients as given from the least square fit. 


A. Pure gravity sector 

First, let us focus on the pure gravity sector alone and 
neglect the (d^g)"^ coefficients. It has been shown in [15] 
that the main oscillations in the radial distance between 
the Earth and the Moon due to the coefficients de¬ 
pend on 6 linear combinations: 

500 c and 500 s- They can be expressed in terms of the 
standard SME coefficients expressed in an Earth equato¬ 
rial frame and in terms of the longitude of the ascending 
node a and of the inclination P of the Moon’s orbit with 
respect to the equator. These combinations are given by 
Eqs. (107-108) from [15]. The longitude of the ascending 
node a with respect to the equator oscillates around 0. 
This oscillation is due to the secular advance of the lon¬ 
gitude of the ascending node with respect to the ecliptic. 
Similarly, the inclination of the Moon’s orbit with re¬ 
spect to the equator oscillates around /? = 23.44 °. As a 
consequence, the transformation of the LLR linear com¬ 
binations to the standard SME coefficients is given by 


■Sllr 

*LLR 

*LLR 

5lLR 

5lLR 

*LLR 


gll _ ^ 0 92 

(16a) 

+0.08 - 2s^^) - 

0.73s^^, 

5^2 = 0.92s^^ + 0.40s^^ , 

(16b) 

= 0.92s^’^ + 0.40s^^ , 

(16c) 

, 

(16d) 

sn^c = -3.21s^^ - 1.39s^^ , 

(16e) 

sq^s = -3.50s^^ . 

(16f) 


Note that the above transformations are different from 
those used in [21]. In that paper, the authors have used 
a = 125 °, which corresponds to the transformation be¬ 
tween the lunar plane and the ecliptic plane at the date 
J2000 while the reference frame used in the SME frame¬ 
work is the equatorial plane (and not the ecliptic one). 
Therefore, the value of a and /3 needs to be taken with 
respect to the equatorial plane at the moment where the 
experiment was performed, or as their average value if 
they vary during the experiment.^ 

In [19], Battat et al have fitted the amplitudes re¬ 
lated to the signature of the 6 SME combinations (16) 
on residuals of LLR analysis. As a result, they obtained 
constraints given in Tab. V. 


^ Note that [15] advised caution on this point: “For definite¬ 
ness and to acquire insight, we adopt the values a = 125 ” and 
p = 23.5 °. However, these angles vary for the Moon due to 
comparatively large Newtonian perturbations, so some caution 
is needed in using the equations that follow.” 
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TABLE V. Estimations of the SME coefficients derived from 
LLR analysis from [19]. 


SME linear combination 

Estimation 

Sllr 

(1.3 ±0.9) X 10-1° 

Sllr 

(6.9 ±4.5) X 10-11 

Sllr 

(-5.2 ±4.8) X 10-1’ 

Sllr 

(-0.8 ± 1.1) X 10-° 

Sllr 

(0.2 ±3.9) X lO-i 

Sllr 

(-1.3 ±4.1) X lO-i" 


B. Gravity sector and matter-gravity couplings 


Combining these constraints with those obtained in the 
previous section from planetary ephemerides lead to es¬ 
timations of the pure gravity SME coefficients given in 
Tab. VI. One can see that the and the three 

coefficients (with J ^ K) are improved by the com¬ 
binations of the data. This is mainly due to the fact that 
the correlations are reduced. It is also worth mentioning 
that this combined analysis improves the combined LLR 
and atom interferometry gravimetry analysis from [21] 
by 2 to 3 orders of magnitude. 


TABLE VI. Estimated mean and la uncertainty of the SME 
coefficients by combining planetary ephemerides analysis 
from Sec. Ill and LLR analysis [19]. It has been assumed that 
the coefficients vanish. 


SME coefficients 

Estimation 

- si'i" 

(9.6 ±5.6) X IQ-ii 

gO = - 2 

(1.6 ±0.78) X 10-1° 


(6.5 ±3.2) X IQ-ii 


(2.0 ± 1.0) X 10-11 


(4.1 ± 5.0) X 10-1° 

-TX 

(4.3 ±2.5) X 10-° 

^TY 

(1.1 ± 1.1) X 10-® 

-TZ 

(-3.8 ±3.0) X 10-® 


In order to use LLR analysis to constrain simultane¬ 
ously the and (a^g)'^ coefficients, we need to identify 
the contributions of the coefficients to the ampli¬ 

tudes of the Earth-Moon distance oscillations. The SME 
contribution to the equations of motion of the Moon- 
Earth system can be found in [15, 18] and is given by 


2^J 


d^r 


dt'^ 


GnM 


SME 


^JK^k _ 


L j,K j,L rpJ 


-f 3 


-TK_ 


E 

w—e,p,n 


5 f « (a-,)'' ) - V’fXrX (17) 




^6m 


-TK 




TV 1 

^ - 2^ 
M 


7:TJ 


E 

w—e,p,n 




where M = m(j -I- TOq, Sm = TOq — TO(j « TOq, is 
the position of the Moon with respect to the Earth, 
is the relative velocity of the Moon with respect to the 
Earth, is the heliocentric velocity of the Earth-Moon 


Barycenter and 


= Nl 


s - Vi" 

" -''e > 

(18a) 

1 




= M 

— > 

(18b) 

\ 


"m-®) 


where is the number of particles of species w in the 
body i. Following the approach described in Appendix 
A of [15] (see also [46, 47]), we expand the equations of 



















motion around a reference circular orbit and perform a 
Fourier analysis to obtain the contributions of the 
terms to the oscillations of the Earth-Moon distance. The 
term proportional to in the first line of Eq. (17) 

leads to an oscillation at the Earth orbital frequency . 
The (dgfi) coefficient modifies the expression of 
in Eq. (A20) from [15]. Similarly, the modifications of 
the terms proportional to Sm in Eq. (17) change the ex¬ 
pression for and To summarize, we find that 

the (a^fg)"^ coefficients will modify the combinations ap¬ 
pearing in LLR oscillations as (s^li^ and being un¬ 
changed) 


Sllr - + X! 


_ /v^ 

-02 ^ /-w \2 


(19a) 


Sllr = 


E '^2 ^ /-w 
w 

W ^ 

^LLR = + 2 ^ ^cos 77 a (a“g) -I- sin r]a (a“g) 

W 

+ + (cosr7a«ff)^ 


(19b) 


= sn„ 


-hsin 77 a(a“ff) ) ,(19c) 


*LLR 


UJ 

= snes + 2^^a«ff)^ 


+ S (iM) 


where sq^c and sn^s are given by Eq. (108) of [15] or 
by Eqs. (16e-16f). A simple model for the composition 
of the Earth leads to IV^/toq = N^/m^ « N^/m^ « 
0.5(GeV/c^) ^ [18]. Similarly, the model for the com¬ 
position of the Moon from [48] leads to A(| /m(j = 

/to(j ~N^/m(^ «0.5(GeV/c^) Using these val¬ 
ues, the combinations (16c-16f) appearing in LLR data 
analysis are modified by the coefficients as follow: 


s^:lr = 0.92 - 0.5a (a:+^)'' - 0.5a (a^g)"’) (20a) 

-1-0.4 - 0.5a (oeff^)^ - 0.5a (dgg)^'j , 

Sllr = - 0-5a {dl'^^)^ - 0.5a (a”g)^ , (20b) 

s^lr = -3.21s^'^ - 1.39s^^ -F 1.84a (20c) 

-|-1.84a (a"g) -I-0.8a (og^^) -I-0.8a ( 0 "^) , 

sLk = -3.50s-^^ + 2a {dl^^f + 2a (hg^)^ . (20d) 

Atom interferometry gravimetry has also been used 
to constrain SME coefficients [20, 21]. A violation of 


Lorentz symmetry induces periodic variations of the lo¬ 
cal acceleration that can be measured by atom gravime¬ 
try. Amplitudes of these oscillations have been partially 
computed in [15] for the coefficients (see Table IV) 
and in [18] for the (d^g)'^ coefficients (see Table IV). An 
improved calculation shows that the (d^g)"^ coefficients 
modify only two of the amplitudes constrained in [ 20 , 21 ]: 


a = *i,-^^sin2x-2^75a-^^ 
2 c 


4Vl 


w—e,p,n 


3V" 




rriT 2 77i0 




*4 -A • o 

= YSAiSm2x, 

*4 7;TX 


( 21 a) 


Dui = —s sin2x-|-2 —755 

1 c 


3c ^ 


m 


3V11 




w—e,p,n 


rriT 2 




X 


^4 • o 

= ^SAiSm2x, 


( 21 b) 


where % = ~ 1/2 (with the Earth 

spherical inertial moment and the Earth radius), 
*4 = 1 —3i0 « —1/2, ib = 1+270/3 « 4/3, the subscripts 
T refer to the test body, Vl = siny is the velocity 

of the laboratory due to Earth rotation (a ;0 being the 
angular velocity of the Earth rotation) and x is the geo¬ 
graphical colatitude of the location where the experiment 
is performed. In the last expressions, we introduced two 
linear combinations given by 


4 


+ 


74 sin 2 x c 

8Vl 1 

3c 74 sin 2x 


Ul . -TY 
75 s 


E 


sfl = + T 


w—e,p,n 

4 Vl. 


74 sin 2 x c 
8Ul 1 




E 


3c 74 sin 2x ^ 

w—e,p,rL 





( 22 a) 


+ 3 /^' 

a 


rriT 

2 7770 _ 





( 22 b) 


+ 3 ^' 

a 

■ 

rriT 

2 7770 



For the experiment performed by [20, 21], we have x = 
42.3 ° and Vl/c « 1.04 x 10“®. Moreover, numerical 
estimations for a Gaesium atom interferometer lead to 
^Cs/"^Cs = N^Jmcs = 0 . 44 (GeV/c 2 )“\ Ngjmcs = 
0.63(GeV/c^) Finally, the values for the Earth are 
given in [18] and are mentioned above after Eq. (19). 
Using these values gives 


-p 1.12 X 10-^s^’^ (23a) 

-5.43 X 10-®a (cigj^)^ - 5.96 x lO'S {d'^gf , 
sfj = - 1.12 X 10-^s^^ (23b) 

-P5.43 X 10-®a (cigff^)^ + 5.96 x lO'^a (a”g)^ , 


with (ogg^)'^ given by Eq. (6). 
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Therefore, the experiment from [20, 21] is sensitive to 
the last two combinations and not to and alone. 
The results from [21] are presented in Tab. Vll 

TABLE VII. Estimations of the SME coefficients derived from 
atom interferometry gravimetry by [20, 21]. 


SME linear combination 

Estimation 


(4.4 ± 11) X 10"® 


(0.2 ±3.9) X 10"® 

-5aI 

(-2.6 ±4.4) X 10"® 

7;B 

Sai 

(-0.3 ±4.5) X 10"^ 

-TX 

(-3.1 ±5.1) X 10“® 


(0.1 ± 5.4) X 10"® 

^TZ 

(1.4 ± 6.6) X 10"® 


In our final analysis, we combine the three analysis 
with both the and coefficients: (i) plane¬ 

tary ephemerides analysis given by Tab. II with the cor¬ 
relation matrix from Tab. Ill (or equivalently the re¬ 
sults from Tab. II on the linear combinations given by 
Eqs. (15)), (ii) LLR data analysis from [19] summarized 
in Tab. V with linear combinations given by Eqs. (16a- 
16b) and (20) and (iii) atom interferometry gravimetry 
analysis from [20, 21] presented in Tab. VII with the lin¬ 
ear combinations given by Eq. (23). The (marginalized) 
results of this fit are presented in Tab. VIII. 

TABLE VIII. Estimated mean and la uncertainty of the SME 
coefficients obtained with a fit combining results from Sec. Ill, 
LLR data analysis from [19] and atom interferometry gravime¬ 
try experiment [20, 21]. 


SME coefficients 

Estimation 


-XX _ -VV 

(9.6 ±5.6) X 10"“ 


sQ = ± - 2 

(1.6 ±0.78) X 10"^° 



(6.5 ±3.2) X 10"“ 



(2.0 ± 1.0) X 10"“ 



(4.1 ± 5.0) X 10"^^ 


-TX 

(-7.4 ±8.7) X 10"® 



(-0.8 ±2.5) X 10"® 


-TZ 

(0.8 ±5.8) X 10"® 



(-7.6 ±9.0) X 10"® 

GeV/c® 


(-6.2 ±9.5) X 10"® 

GeV/c® 

Q! (ffllfl)^ ± O (“eff)^ 

(1.3 ±2.2) X 10"'‘ 

GeV/c® 

a(aSr)"' 

(-5.4 ±6.3) X 10"® 

GeV/c® 

«(a:ff)^ 

(4.8 ±8.2) X 10"'‘ 

GeV/c® 


(-1.1 ± 1.9) X 10"® 

GeV/c® 


The resulting estimations do not show any significant 
deviations from GR. The combinations of the three data 
analyses allow one to estimate each of the coefficients 
individually. The spatial part of is completely de¬ 
termined by the combination of planetary ephemerides 
and LLR data. The atom interferometry gravimetry is 


not accurate enough to provide any significative improve¬ 
ment on the uncertainty of these coefficients. With an 
improvement of 2 orders of magnitude, the atom gravime¬ 
try data would become significative to estimate the 
coefficients. On the other hand, the three datasets are re¬ 
quired in order to decorrelate the s'^'^ and the (d[][j)'^ co¬ 
efficients. The uncertainties on are much larger than 
those shown in Tab. VI where the coefficients (d((ff)'^ have 
been neglected. This reflects the fact that the individual 
coefficients are still highly correlated. 


V. DISCUSSION 

First of all, the accuracy of the constraints on the SME 
coefficients obtained in Tab. II (planetary orbital dynam¬ 
ics alone) are of the same order of magnitude as the bi¬ 
nary pulsars [25] constraints on the SME coefficients with 
an improvement of one order of magnitude on the coeffi¬ 
cients . Nevertheless, it is known that nonperturba- 
tive effects (similar to those computed in [49]) may arise 
in binary pulsar systems. The nonperturbative effects de¬ 
pend highly on the fundamental theory (for example, see 
[50] for nonpertubative calculations in Einstein-Aether 
theory or in Hofava gravity). In general, the results from 
[25] are effective constraints on the strong field version 
of the that may include nonperturbative strong field 
effects and one should be careful when comparing strong 
field tests and weak field tests as the one performed in 
Sec. III. The results shown in Tab. IV improve the cur¬ 
rent Solar System constraints [32] by I to 3 orders of 
magnitude. Furthermore, the analysis combining plane¬ 
tary orbital dynamics and LLR from Tab. VI improves 
by 2 to 3 orders of magnitude the previous results that 
combined LLR and atom interferometry. This shows the 
high impact provided by planetary ephemerides analysis. 

As mentioned in Sec. Ill, our results show that the 
estimated SME coefficients are highly correlated. The 
correlations are due to the similarity of the orbital planes 
of all the planets. Therefore, one way to improve the 
results by reducing the correlations is to use bodies with 
different orbital planes like e.g. asteroids. This can be 
achieved for example with Gaia observations similar to 
what is proposed in [51]. 

The constraints obtained in Sec. Ill are mainly due to 
the internal planets. For instance, Jupiter has absolutely 
no influence on the results shown in Tab. II. This is a 
consequence of its not so well-known orbit. An improve¬ 
ment by a factor 10 on the knowledge of Jupiter’s orbit is 
required for that planet to play a significant role in this 
analysis. Therefore, the improvement of Jupiter’s tra¬ 
jectory expected from the analysis of Juno’s radioscience 
and very long baseline interferometry data [52] may im¬ 
prove the result of our analysis. In particular, it will 
reduce some of the correlations which will lead to an im¬ 
provement of the estimations of the SME coefficients. In 
the same spirit, the influence of Saturn is weak but nev¬ 
ertheless highly important to decorrelate the coefficients. 
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Furthermore, an improvement of Mercury’s orbit by a 
factor 10 (which can be regarded as the improvement by 
Messenger’s data that are not yet included in INPOPlOa 
analysis [33]) will lead to an improvement on the estima¬ 
tions of by a factor 2 and to a 10 % improvement on 
the coefficients S'^ and (but to no improvement at 
all on the other coefficients). In summary, the best way 
to improve the current analysis is to improve the trajec¬ 
tory of the “badly” determined planetary orbits in order 
to improve the decorrelation instead of improving more 
the planets that are already very well determined. 

As mentioned in Sec. II, the influence of the and 
the coefficients on the orbital dynamics only ap¬ 

pears at the next post-Newtonian order and these coeffi¬ 
cients are therefore not constrained by our analysis. Nev¬ 
ertheless, these coefficients will play an important role 
in the light propagation [53, 54]. Therefore, planetary 
ephemerides may potentially constrain this coefficient by 
considering the effect of on the light time of the ra¬ 
dioscience Range observables used in the analysis. Other 
opportunities to constrain this coefficient are to consider 
a conjunction experiment like the one performed with the 
Cassini spacecraft [3] (or to analyze Cassini data within 
the SME formalism as proposed in [55]) or to consider 
Very Long Baseline Interferometry observations similar 
to what has been done for the 7 post-Newtonian param¬ 
eter [5]. 

The multiplication of the numbers of SME coefficients 
that need to be considered leads to an increase in the 
uncertainties on each individual coefficients. This is due 
to the correlations between the different coefficients that 
appear when their number is increased. Therefore, it is 
highly important to increase the number of analyses to 
constrain SME. In this communication, we have shown 
how a combination of three analyses can disentangle the 
different coefficients. Nevertheless, the coefficients shown 
in Tab. VIII are still highly correlated, especially in the 
(a^)"^ sector. One way to reduce these correlations is to 
use more observations that are sensitive to other combi¬ 
nations of the {a^ffY coefficients. This can be done in two 
ways: (i) to consider different source bodies that generate 
the gravitational field and (ii) to use more orbital geom¬ 
etry like e.g. asteroids dynamics as already mentioned. 
The first point is related to the fact that the (a]([j)'^ coef¬ 
ficients enter the equations of motion essentially through 
the properties of the source body. In this communica¬ 
tion, only two source bodies have been used: the Sun (in 
the planetary orbital dynamics analysis) and the Earth 
(in LLR and in atom interferometry gravimetry). Con¬ 
sidering more source bodies with different compositions 
can help to reduce correlations. In this sense, a test us¬ 
ing the satellites around the different planets would be 
highly relevant. 

Einally, we would like to soften the results presented 
here. First of all, we insist on the fact that the con¬ 
straints obtained in Sec. Ill correspond to the intervals 
in which the differences of INPOPlOa postfit residuals are 
below 5 %, as they are obtained directly from the limits 


of Tab. I coming from [33]. As such, they do not directly 
represent the usual Icr confidence interval. A cleaner ap¬ 
proach would be to include the SME equations of motion 
directly in the planetary ephemerides software and to es¬ 
timate the SME coefficients directly from the raw data, 
which corresponds to the approach usually used for esti¬ 
mating the PPN coefficients [4, 6 , 7] or more recently to 
constrain the MOND theory [56]. Our analysis demon¬ 
strates the impact of such an analysis and therefore, pro¬ 
vides a strong incentive. 

In addition, the LLR data analysis has been performed 
by fitting some oscillating signatures in the LLR data 
residuals. This approach is not optimal since it suffers 
from two drawbacks. First, the oscillating signatures de¬ 
rived in [15] have been computed analytically using sev¬ 
eral approximations. They can be used to estimate an 
order of magnitude on the different effects produced by 
SME but they are not optimal for a real data analysis 
(furthermore, the signatures used in [15] include only the 
dominant oscillations, several other frequencies are pro¬ 
duced by SME and ignored in the data analysis). Second, 
fitting in the residuals is not optimal since it does not 
allow one to analyze the correlations between the SME 
coefficients and the other parameters that are usually fit¬ 
ted in a standard LLR data analysis. Eor these reasons, 
a cleaner analysis would include the SME equations of 
motion directly in the software used to reduce LLR data. 
Results obtained in [19] and in this communication give 
strong motivations to perform such an analysis. 

Finally, the atom interferometry gravimetry analysis 
should be interpreted with caution. The atom interfer¬ 
ometry gravimeter results from [ 20 , 21 ] assume a model 
of the local solid Earth tides. While such models can be 
partly analytically based, it is known that the many fre¬ 
quencies of the Earth tides include all of the frequencies 
in the SME signal [57]. If any aspect of the tidal model 
includes fitting sinusoidal functions to local gravimetry 
measurements or global measurements of, for example, 
the ocean heights [58], the signal for the SME may be 
partly subtracted due to the strong correlation with the 
tidal signal. 


VI. CONCLUSION 

In this communication, we have shown that the plan¬ 
etary orbital dynamics allow one to constrain a viola¬ 
tion of Lorentz symmetry with an impressive accuracy. 
In Sec. Ill, we use the current limits on supplemen¬ 
tary advances of perihelia and nodes provided by IN¬ 
POPlOa [33] to estimate the SME coefficients and 
{a^gY. In this analysis, the coefficients have been 
neglected since they are already constrained with a high 
level of accuracy [32] but they can be considered in a 
future work. Our analysis has been performed using a 
standard Bayesian inversion. Results on the SME coef¬ 
ficients are given in Tab. II. No significative deviation 
from GR is observed. As mentioned in Sec. Ill, these 
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estimations are highly correlated (see Tab. Ill or Fig. 1). 
We have identihed numerically the linear combinations 
of the SME coefficients that can be estimated indepen¬ 
dently from planetary ephemerides. The estimations on 
these combinations are given in Tab. IV. These two re¬ 
sults are completely equivalent (as long as one uses the 
correlation matrix with the first estimation). Our results 
produce uncertainties similar to those obtained from bi¬ 
nary pulsars data [25] on most of the coefficients and 
improve the constraints on by one order of magni¬ 
tude. Moreover, we improve the current best weak field 
tests by 2 to 3 orders of magnitude. 

We also perform a combined estimation of the SME co¬ 
efficients using results from three different analyses: (i) 
the planetary ephemerides analysis performed in Sec. Ill, 
(ii) the LLR data analysis performed in [19] and (hi) 
the atom interferometry gravimetry analysis realized 
in [20, 21]. The combination of LLR and planetary 
ephemerides leads to the best current estimations on 
the pure gravity SME coefficients as shown in Tab. VI 
(when neglecting the (d[(g)‘^ coefficients). In these three 
analyses, we also take into account potential effects pro¬ 
duced by a Lorentz violation in the matter-gravity cou¬ 


pling which is parametrized by the coefficients. 

Einally, the combination of the results from the three 
data analyses leads to the first independent estimations 
of the and (a^g)'^ coefficients. The results are pre¬ 
sented in Tab. VIII. The obtained uncertainties are rel¬ 
atively large, which is due to the numbers of coefficients 
considered and to the remaining correlations. Some ideas 
to reduce these correlations are proposed in Sec. V. 
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